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Abstract 

The pion and kaon electromagnetic form factors Fm{Q'^) are calculated 
at the leading order of pQCD using the running coupling constant method. 
In calculations the leading and next-to- leading order terms in Q5((1 — — 
y)Q'^) expansion in terms of asiQ"^) are taken into account. The resummed 
expression for Fm{Q'^) is found. Results of numerical calculations for the 
pion (asymptotic distribution amplitude) are presented. 
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During last years a considerable progress were made in understanding of infrared 
(ir) renormalon effects in various inclusive and exclusive processes [1-4]. It is well 
known that all-order resummation of ir renormalons corresponds to the calculation 
of the one-loop Feynman diagrams with the running coupling constant as{—k'^) 
at the vertices, or to calculation of the same diagrams with non-zero gluon mass. 
Both these approaches are generalization of the Brodsky, Lepage and Mackenzie 
(BLM) scale-setting method [5]. Studies of ir renormalon effects have also opened 
interesting prospects for evalution of higher twist corrections to processes' different 
characterictics [6]. 

As was proved in our works [2-4], exclusive processes have additional source of ir 
renormalon corrections. Indeed, integration over longitudinal fractional momenta 
of hadron constituents in the expression, for example, of the electromagnetic (elm) 
form factor generates ir renormalon effects. 

In this letter we extend our previous consideration of the pion and kaon elm 
form factors in the context of the running coupling constant method. 

A meson M form factor in the framework of pQCD has the following form [7] 
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dxdy(j)*M (y, Q^) Th (x, y; Q^, as (Q^)) 0m {x, Q^) . (1) 



Here — —q^ is the square of the virtual photon's four-momentum, (pM {x, Q^) 
is the meson distribution amplitude. In Eq.(l) Th (x,y;Q'^,as is the hard- 

scattering amplitude of the subprocess qq' -|- 7* — > qq', which at the leading order 
of pQCD is [7] 
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In Eq.(2) Cp = | is the color factor, is taken as the square of the momen- 
tum transfer of the exchanged hard gluon in corresponding Feynman diagrams 
for Fm (Q^) ■ Unlike the frozen coupling approximation, where for calculation of 
Fm (Q^) one fixes the argument of 0:5 ((1 — x) (1 — y)Q'^) as (Q^/4) , in the run- 
ning coupling method we use Eq.(2) in Eq.(l) and, as a result, encounter with 
new troubles. Indeed, it is evident that as (Q^) Eq.(2) suffers from ir singular- 
ities associated with the behaviour of the as (Q^) the regions x ^ 0,y — 0; 
a; ^ 1,1/ — > 1. Therefore, Fm {Q"^) can be found only after proper regularization of 
o^s (Q^) ill these end-point regions. For solving of this problem let us relate the run- 
ning coupling as (AQ^) in terms of as (Q^) by means of the renormalization group 
equation. The renormalization group equation for the running coupling a = as/'K 
has the form 
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where 



&2 = ^ (33 - 2n;) , 63 = ^ (306 - 38nf) . 
The solution of Eq.(3) with initial condition a{X) \x=i— a = as{Q'^)/7r is [8] 
a{X) 
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This transcendental equation can be solved iteratively by keeping the leading In '^A 
and next-to- leading In *^~^A powers. For A = [1 — x){l — y) these terms are given 

by 



-(^^-^^^-^^^^)- i + .(g J- .)(!-.) 

Q;^(Q^)b3 In [1 + a{Q'')h2 ln(l - x){l - y)] 
62[l + a(g2)62ln(l-x)(l-y)]' 

The first term in Eq.(5) is the solution of the renormalization group equation (3) 
with leading power accuracy, whereas the whole expression (5) is the solution of 
Eq.(3) with next-to- leading power accuracy. 

In our previous papers [3,4] for calculation of the pion and kaon elm form factors 
we used only the first term from Eq.(5). In this work for evaluation of FuiQ"^) we 
use both of them. Let us clarify our approach by considering the pion form factor 
and pion's simplest distribution amplitude (f)asy{x) 

(Pasy{x) = V3Ux{l - X). (6) 

Generalization of obtained results for the pion's other distribution amplitudes as 
well as for the kaon is straightforward. 

After substitution of Eqs.(5),(6) into Eq.(l), the pion form factor takes the form 
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where t — l/ab2- 

The resummed expression for Ft^{Q'^) can be found using an approach advocated 
in Ref. [3]. For these purposes it is instructive to change variables from y in Eq.(7) 
to ^ = ln(l — x),w = ln(l — y) and apply the inverse Laplace transforms [9] 
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and 



hi (t + z + w) 
(t + z + wf 



— / exp[—u(t + z + w)](l — C — ln.u)udu. 
Jo 
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In Eq.(9) C ~ 0.577216 is the Euler-Mascheroni constant. After some manipula- 
tions one gets 
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exp{-tu)B Q^Ft, {u)du 



ex^{-tu)u[l-C -hit-\nu\B Q^F^ {u)du, (10) 



where B [Q^F^] {u) is given by expression 
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and is the Borel transform of the pion elm form factor obtained in our work [3]. 
It is worth noting that Eq.(lO) is the general expression valid for both the pion 
and kaon; one needs only to replace in Eq.(lO) the Borel transform B [Q'^F^^] (u) 
and with corresponding ones B [Q'^F^] (u) and fx from Ref.[4]. The first term 
in Eq.(lO) is the pion form factor found in Ref.[3] at the leading order of pQCD 
in the framework of the running coupling method using only the first term from 
Eq.(5). Of course, the second term in Eq.(lO) is the contribution to the form factor 
coming from the second term of Eq.(5). But interpretation of the integrand in 
this term (without exp(— tw)) as a "traditional" Borel transform of corresponding 
perturbative series is problematical, because some regularization prescription at 
M = for recovering of perturbative series is needed. Here we bypass this problem 
and concentrate on calculation of the resummed expression for Fm{Q^)- As we shall 
see later, integration in Eq.(lO) using the principal value prescription [8,10] removes 
ir renormalon poles at u— 1, 2... and gives correct results for [Q'^Fm{Q'^)Y^^ ■ 
For the pion we get 
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Here [Q'^ F-,^{Q'^)\'^^ is the resummed elm form factor of the pion obtained in Ref.[3] 
by means of the leading term of Eq.(5), 
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In Eqs.(12,13) we introduce following notations, 
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For the pion's Chernyak-Zhitnitsky (CZ) distribution amplitude [11] in the form 
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we get 
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Expressions for the resummed pion form factor obtained using more general than 
the asymptotic and CZ distribution amplitudes [12] are rather lengthy and will be 
published elsewhere. 
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For the kaon we find 



+ (C + lni)mfc/fc(0) - {tm.k + nk)gk{l) + 111^5(^(0)} . 

The values of (rrifc, n^) as well as [Q^ Fx{Q'^)\'i^ can be found in Ref.[4]. 

As an example, results of numerical calculations for [Q"^ Ft^{Q'^)Y'^^ carried out 
using the pion's asymptotic distribution amplitude are shown in Fig.l. Here the 
ratio R = [Q^F,r(Q^)]'''V[<5^i^^(<5^)]°, where is the pion elm form 

factor calculated in the frozen coupling approximation using the same distribution 
amplitude, is depicted. As is seen, the next-to-leading order correction changes the 
shape of the curve in the region of small (2 GeV^ < < 6 GeV"^), where the 
whole result is smaller than the leading order one. At = 2 GeV'^ the correction 
amounts to ~ 15% of the leading order result decreasing toward Q"^ — 6 GeV"^. 
At the end of the considered values of the next-to-leading order correction is 
positive and equals to 3 — 6% of the leading order result. The similar calculations 
can be also fulfilled for the pion using its alternative distribution amplitudes [12] 
and for the kaon. 

As was pointed out in Refs.[3-5], the ir renormalon corrections can be hidden 
into the scale of asiQ"^) at the leading order elm form factor [Q^Fm{Q^)Y'- In the 
studied case of the pion (asymptotic distribution amplitude) we find 

f{Q^)c^-6.71 + 16.67as{Q^), (15) 

where in numerical fitting we have used Eq.(12). 

It is known [3,4] that ir renormalon effects enhance the perturbative predictions 
for the pion, kaon elm form factors approximately two times. Our recent stud- 
ies confirm that next-to- leading order term in Eq.(5) does not change the picture 
considerably. This feature of ir renormalons may help one in solution of a contradic- 
tion between theoretical interpretations of experimental data for the photon-to-pion 
transition form factor F^^^^Q'^) [13] from one side and for the pion elm form fac- 
tor Ft^{Q'^) [14] from another side. Thus, in works [15,16] the authors noted that 
the scaling and normalization of the photon-to-pion transition form factor tends 
to favor of the pion asymptotic-like distribution amplitude. But then prediction 
for Ft^{Q^) obtained using the same distribution amplitude is lower than the data 
by approximately a factor of 2. We think that in this discussion a crucial point 
is a chosen method of integration in Eq.(l). Indeed, unlike F^^IQ"^) the expression 
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for F^7r(<5^) at the leading order of pQCD does not contain an integration over 
asixQ"^). In other words, the running couphng method being apphed to Eq.(l) and 
to the expression for (see, Refs. [15,16]) enhances the perturbative result for the 
pion elm form factor and, at the same time, does not change -F^jrQ- This allows us 
to suppose that in the context of pQCD the same pion distribution amplitude may 
explain experimental data for both F^t^^Q"^) and Ft^{Q'^). This problem is a subject 
of separate publication. 



^Our statements concerning F^^{Q'^) remain true also in the light of ir renormalon corrections 
computed for F^.^ in Ref. [17], because they have another source than ones considered in our work, 
namely, Feynman diagrams with bubble chains. At the leading order of pQCD the only source 
of ir renormalons is the running coupling constant and integration in Eq.(l) over longitudinal 
momenta of a meson M . 
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FIGURE CAPTION 

Fig.l The ratio R = [Q^F^{Q^)Y^' /[Q^F^{Q^)]^ as a function of Q^. in calculation 
the pion asymptotic distribution amplitude (pasyi^) is used. The curve 1 is i? found 
using the whole result for the resummed form factor Eq.(12), whereas the dashed 
curve is the ratio R — [(3^-^7r(Q^)]rV['5^-^7r(Q^)]°- The curve 2 corresponds to 
R with [Q^ Ft^{Q'^)]'^ in the frozen couphng approximation and after scale-setting 
procedure (15). In calculations the QCD scale parameter A has been taken equal 
to 0.1 GeV. 
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